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Abstract: Entanglement entropy of holographic CFTs is expected to play a crucial role 
in the reconstruction of semiclassical bulk gravity. We consider the entanglement entropy 
of spherical regions of vacuum, which is known to contain universal contributions. After 
perturbing the GET with a relevant scalar operator, also the first order change of this quantity 
gives a universal term which only depends on a discrete set of basic GET parameters. We 
show that in gravity this statement corresponds to the uniqueness of the ghost-free graviton 
propagator on an AdS background geometry. While the gravitational dynamics in this context 
contains little information about the structure of the bulk theory, there is a discrete set of 
dimensionless parameters of the theory which determines the entanglement entropy. We argue 
that for every (not necessarily holographic) GET, any reasonable gravity model can be used to 
compute this particular entanglement entropy. We elucidate how this statement is consistent 
with AdS/GFT and also give various generalizations. On the one hand this illustrates the 
remarkable usefulness of geometric concepts for understanding entanglement in general GFTs. 
On the other hand, it provides hints as to what entanglement data can be expected to provide 
enough information to distinguish, e.g., bulk theories with different higher curvature couplings. 
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1 Introduction 

It is by now well known that some conformal field theories (CFTs) admit a dual description 
in terms of semiclassical gravity in one higher dimension [1]. But it is so far not clear what 
are the precise CFT data that one needs to know in order to reconstruct various specific 
features of the bulk theory. A full knowledge of the CFT should be exactly equivalent to full 
knowledge of bulk quantum gravity, the challenge being to find the precise dictionary. We 
wish to focus on the question which subset of CFT parameters one needs to have access to 
in order to reconstruct just a classical bulk geometry and its linearized dynamics. Once this 
question is answered in terms of fundamental CFT parameters, the practical question is which 
CFT observables are sensitive to these parameters in a useful way. 
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One hint is the fact that a local bulk metric should be represented in the boundary CFT 
in a non-local fashion. An important CFT feature which does have non-local properties is 
the entangled structure of the underlying quantum state. It has therefore been proposed that 
knowing enough about the entanglement structure of the CFT state is sufficient to determine 
the gravitational dynamics of the bulk. In particular, entanglement entropy of spatial regions 
in the CFT appears to be an illuminating quantity to study in this context, c.f., [2-6]. The 
basic rationale for this idea is that entanglement entropy of a given spatial region A can be 
computed in semiclassical holography as the area of a certain extremal surface [7-9] which, 
of course, encodes information about the metric. By varying this entanglement entropy (for 
instance by changing A, or by varying the underlying quantum state, or by deforming the 
CFT spectrum), the corresponding change of the minimal surface encodes information about 
the dynamics of the bulk metric. Making this idea precise and general would be significant 
progress towards our understanding of the holographic nature of quantum states and gravity. 

We wish to carefully distinguish differently strong versions of the expectation that en¬ 
tanglement encodes bulk geometry. The weakest version of the “entanglement is geometry” 
statement simply says that the most efficient way of calculating entanglement entropy is by 
using geometric methods. For holographic CFTs the paradigmatic manifestation of this state¬ 
ment is the Ryu-Takayanagi conjecture [7, 8]. But even for non-holographic CFTs there is 
evidence that geometric concepts sometimes provide the most natural and efficient way of 
computing entanglement entropy. The first such statement is the realization that vacuum en¬ 
tanglement entropy of spherical regions can be conformally mapped to thermal entropy which 
sometimes has a natural interpretation in terms of black hole thermodynamics [10]. Refining 
this idea, Faulkner [11] revisited the problem of computing the response of vacuum entangle¬ 
ment entropy of spheres to a deformation of the CFT by a relevant scalar operator [12-17]. 
For this computation, nothing is assumed about the existence or the nature of a holographic 
dual of the CFT. As a CFT calculation this is a formidable problem. Intriguingly, the most 
convenient way of repackaging this calculation turns out to be in terms of an auxiliary grav¬ 
itational problem, where the entanglement entropy is represented as a minimal surface area 
in cosmological Einstein gravity responding to a scalar field perturbation. This is a power¬ 
ful result, in particular in face of the fact that entanglement entropy is notoriously hard to 
compute otherwise. One of the goals of this paper is to elucidate from a gravitational point 
of view how this universality comes about and to reconcile these statements with AdS/CFT 
duality (where bulk descriptions generically have a structure which is not necessarily given by 
just Einstein gravity). 

A stronger form of the “entanglement is geometry” conjecture is the idea that in holo¬ 
graphic CETs the detailed bulk dynamics can be extracted entirely from knowledge of entan¬ 
glement entropy for a suitable set of regions A. Even in the semiclassical regime, this is a 
strong assertion because the detailed bulk dynamics may involve not only Einstein gravity, 
but also higher curvature corrections [18-24]. The imprint of higher curvature corrections on 
equations of motion should then be reconstructible from CFT entanglement entropies. See 
[4, 6, 25] for recent progress towards deriving linearized Einstein equations from a first law of 
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entanglement [5, 26-28]. It will be a topic of this paper to illuminate the way in which CFT 
entanglement does or does not discriminate between different higher curvature interactions in 
the bulk. 

Let us compare these two versions of the “entanglement is geometry” statement. If spheri¬ 
cal region entanglement entropy in deformations of any arbitrary CFT can be computed using 
linearized Einstein gravity [II], it follows that this quantity is certainly not sufficient to dis¬ 
tinguish even CFTs with a semiclassical gravity dual from those without; let alone distinguish 
a bulk theory governed by Einstein gravity from any other higher derivative theory with the 
same spectrum of low energy excitations. Therefore, for computing this particular type of 
entanglement entropy, geometry is always an excellent tool; but in order to make sure that 
this geometry is actually the lowest order manifestation of a semiclassical AdS/CFT duality, 
more input is required. Of course, the bulk reconstruction from CFT data proceeds gradually: 
the more CFT parameters the entanglement entropies considered are sensitive to, the more 
details of the bulk dynamics are a-priori expected to be derivable from it. 

In this note we wish to investigate how to find entanglement entropies which one can 
expect to be sufficient to reconstruct various bulk features such as linearized dynamics, higher 
curvature couplings etc. In the context of spherical region entanglement entropy and scalar 
operator deformations, we will give a detailed analysis both from a CFT point of view and 
from the bulk perspective. As we will see, at the first relevant order in the perturbation, the 
entanglement entropy is only sensitive to few correlation functions which have a universal form 
(up to normalization) thanks to conformal invariance. From the bulk perspective the analogous 
statement is that the backreaction of a maximally symmetric geometry to the presence of a 
scalar probe is universal for any acceptable theory of gravity.^ A careful comparison reveals 
that this entanglement data is not sufficient to conclude anything about the presence of higher 
curvature terms in the bulk Lagrangian. For the purpose of bulk reconstruction we will outline 
how to gain sensitivity to such features, but will also demonstrate explicitly that the existence 
of some higher curvature interactions is much easier to conclude than the precise structure of 
the latter. 

The structure of this note is as follows. In §2 we first review the context of our discussion 
and then give CFT and gravity arguments to demonstrate universal properties of entangle¬ 
ment entropy of spherically symmetric regions. In §3 we illustrate our abstract discussion 
with the example of curvature squared theories of gravity and Gauss-Bonnet theory in par¬ 
ticular. Finally, §4 contains a discussion and relates to recent developments in this context. 
We generalize our discussion to a wider class of higher curvature theories in §A, allowing also 
for fourth order equations of motion. In §B we perform the most general background field 
expansion for gravity at second order in the metric perturbation, thus uncovering the univer¬ 
sality of Gauss-Bonnet theory within this class of theories. Some facts about the first law of 
entanglement are reviewed in §C. Various useful formulae can be found in §D. 


^ By an “acceptable theory of gravity” we mean any semiclassical gravitational theory with the correct 
spectrum of excitations (in particular without ghosts). See however §A. 
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2 Universality of entanglement entropy of spherical regions 

In this section we analyze universal features both of CFTs and of gravitational theories, which 
explain universal properties of entanglement entropy of ball shaped spatial regions in CFTs 
deformed by a scalar operator. We will then compare these using holography and draw 
conclusions for bulk reconstruction from entanglement. 

2.1 Setup 

Consider a holographic CFT,^ in the vacuum state and on Minkowski spacetime. Assume 
further that the theory has a classical gravity dual described by Einstein theory and consider 
a ball shaped spatial region A of radius R. The reduced density matrix associated with A is 
obtained by tracing out the degrees of freedom outside of A: 

/o.4 = Tiu^|0)(0|. (2.1) 

Without loss of generality, we will assume p_A is normalized, i.e., Tr/9_4 = 1. The entanglement 
entropy of A is dehned as the von Neumann entropy of and can be computed in the bulk 
as the area of the minimal surface ending on the boundary of the ball [7]: 

sf'f M) = -Tr(p.4 log pa ) = , (2.2) 

where Emin is the corresponding bulk minimal surface and Gj^ denotes Newton’s constant in 
d + 1 dimensions. Since the region A is spherical, its causal development can be mapped by a 
conformal transformation to a direct product of time with a maximally symmetric hyperbolic 
space [10]. This transformation further maps the CFT vacuum state to a thermal state. The 
extremal surface is no longer anchored at the boundary but instead wraps the horizon of 
a hyperbolic black hole. The computation of entanglement entropy of A thus reduces to a 
calculation of the horizon entropy of a hyperbolic black hole at temperature T = l/(27ri?): 

= Sthermamd-l) = , ( 2 . 3 ) 

where is the hyperbolic horizon slice. For notational simplicity, we will from now on 

work with Planck units, i.e., = 8ttG^. 

Now let us ask what changes if the bulk gravitational theory is not just Einstein gravity, 
but also contains some higher curvature corrections. For concreteness, consider a d + 1- 
dimensional bulk theory of the form 

where = Gh.dXdab-, Rabcd-,^ eRabcd-, ■ ■ ■) cncodes higher curvature corrections. The cosmo¬ 
logical length scale ^ is related to the AdS scale -^AdSi e-g-i iii Einstein gravity, i = .^AdS- The 
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thermal entropy (and thus the entanglement entropy) is now given by the Wald entropy [29] 
of the hyperbolic horizon [19, 30, 31]: 


= ^r(d/2)a; 


-27r r 
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where Uah denotes the binormal to the hyperbolic horizon normalized as n^'^Uah = — 2 

and 7 is the induced metric. The constant a*^ depends on the details of the gravitational 
theory. For example, in Einstein gravity Ch.d. = 0 and we hnd by comparison with (2.3) 
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In fact, since the Wald entropy density is constant on maximally symmetric AdS^+i, we can 
immediately read off the general expression for the constant 
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Here, the second step is a simple consequence of the existence of an AdS^+i background 
solution and the maximal symmetry of the geometry (see, e.g., [30] for a detailed derivation).^ 
The important feature of (2.5) is that the entanglement entropy is always proportional to 
the horizon area Vol(EIrf_i) (or equivalently to Yol{£rnin)): irrespective of the details of the 
gravitational theory. In particular, just knowing S'gg(^) does not give any insight as to 
what higher derivative corrections the bulk dual of the given CFT may entail: any higher 
derivative corrections to Einstein gravity will only change the overall normalization but not 
the volume scaling in (2.5). But any such change of due to higher curvature corrections can 
be equivalently interpreted as simply a renormalization of the ratio -^Ads/^P within Einstein 
gravity. The precise renormalization factor is given by the bracket in (2.7). The effect of higher 
derivative terms is hence not visible for its only imprint is to renormalize a dimensionless 
parameter of the theory. 

Note that the above expressions are divergent and the matchings performed should hence 
be done more carefully. For example, the volume factor in (2.3) should be regulated such 
that its leading divergent piece scales with the area of the sphere A. This divergence depends 
on the cut-off, but in a consistent way such that the CFT short distance cut-off (regulating 
the UV divergence of entanglement entropy near the boundary of ^) matches the bulk cut¬ 
off (regulating the divergent area of IHIrf_i) [8]. More importantly, the expression (2.3) also 
contains a universal piece which is constant in even dimensions and logarithmically divergent 

^ The first term in the bracket of (2.7) comes from the Einstein contribution in the action (2.4). In a theory 
C = d(d — l)/£^ + Ch.d. without pure Ricci scalar, this term would disappear. 
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in odd dimensions. For this reason, the above matching conditions strictly speaking refer to 
the universal piece, which takes the form [8, 31] 


0(0) 

‘^EE,univ 
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where 5 is the UV cut-off scale in the CFT. It is only for spherical regions A that our arguments 
above give the correct result without explicitly dealing with these additional subtleties, see 
e.g. [19]. 

From a CFT point of view, is a fundamental parameter. As we have just seen, it 
characterizes the normalization of the universal part of vacuum entanglement entropy. It is 
a remarkable fact about entanglement entropy that in even dimensions coincides with the 
coefficient a of the Euler density in the conformal trace anomaly [30-32]. In odd dimensions, 
is proportional to logZgd, the sphere partition function of the CFT [10]. 

For a detailed reconstruction of bulk dynamics from entanglement, one would like to probe 
the particular structure of the gravitational action. In this note we wish to study deformations 
of the above setup and ask what CFT quantities one can use in order to gain sensitivity to 
features which are peculiar to higher curvature terms in the gravity dual. Clearly the vacuum 
entanglement entropy of balls A will not be enough by the above argument. Consider therefore 
a deformation of the CFT by some relevant scalar operator O, uniformly coupled to the CFT 
with a coupling A. Given the arguments above, it is clear what bulk computation we should 
do in order to learn about the change in See{A): add to the gravitational problem a scalar 
field dual to O with suitable boundary conditions, solve the scalar equation of motion in the 
AdSrf_|_i background and compute the backreaction on the minimal surface area perturbatively 
in A. As we will see, the linearized answer is still universal in a way very similar to how (2.5) 
was universal for the unperturbed problem. Namely, its functional form is completely fixed 
such that it only depends on a few parameters which encode the influence of higher curvature 
couplings. These parameters can again be absorbed in a renormalization of dimensionless 
quantities such that the answer is indistinguishable from what one would have obtained in 
Einstein gravity minimally coupled to a scalar field. This is the main universality statement 
studied in this section: For holographically computing entanglement entropy of spherical re¬ 
gions in CFTs deformed by a relevant scalar operator, any physically acceptable theory of a 
spin-2 and a spin-0 field with appropriately tuned values of the couplings gives the same result 
as cosmological Einstein gravity minimally coupled to a scalar. In particular, this quantity 
does not distinguish between different higher eurvature theories of gravity. In the rest of this 
section, we will explain these statements from the CFT and from the gravity side. 


2.2 CFT explanation for universality 

We start by reviewing from a CFT point of view the ingredients involved in computing en¬ 
tanglement entropy of spherical regions in perturbed states. Our discussion closely follows 
[II, 16, 17]. 
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Consider a CFT on d-dimensional Euclidean flat space and perform the following defor¬ 
mation of the CFT action I^prp by a relevant scalar operator O of dimension A < d: 

IcpT = I^cFp + xjd‘^x Oix). (2.9) 

We can compute the perturbative changes of the entanglement entropy See{-^) by writing 
the reduced density matrix as a path integral. This path integral lives on a manifold which 
is the Euclidean spacetime with a cut along dA. The matrix element corresponding to some 
held configurations <pp and (p- on the two sides of the cut is 

(p{dA-)=(p- 

= J^^0A+)=<j>+ e-^clr U-X J d'^xO{x) + yJJ + • • •) 

4>{dA-)=4>- 

( 2 . 10 ) 


Now define a modular Hamiltonian via PA = ^ We can then use the identity 

See{A) = "Tt^paHa) to compute the change in the entanglement to any desired order in 
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J A=0 


= S^^^piA) - X{OH^^) + y - (O^)) + 


( 2 . 11 ) 


where S^pp{A) is universal up to an overall normalization a*^ as explained in §2.1, and similarly 
is the unperturbed modular Hamiltonian. In the last step we used the fact that pA is 
normalized which implies ) = 0. We can now use that Al is a spherical region, which 
brings about another important simplification: the modular Hamiltonian can be traded for an 
integrated stress tensor, since for spherical Al with radius R centered at the origin, one finds 

Ha = 2e ( 2 . 12 ) 

for Al lying on the x® = 0 slice. In our problem, correlation functions involving the modular 
Hamiltonian can therefore be regarded as correlation functions involving the conformal stress 
tensor. But for stress tensor correlation functions, we can use well-known features of CFTs: 
on general grounds, two-point correlations between a stress tensor and a primary operator are 

^ Operators without argument are understood to be integrated, e.g., O = f dfix 0{x). 
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constrained by conformal invariance and conformal Ward identities [33]. In fact, independent 
of the details of our problem, conformal invariance of a theory in a flat spacetime enforces 
{OT^y) = 0. Therefore the 0(A) perturbation in (2.11) vanishes.'^ 

dT 

Furthermore, the perturbation of T^i, is uniform in O and we have = —6^uO. To 
see this, recall that for a Euclidean flat boundary, we have 

T^u{x) = ^ - 5^j_y\0{x). (2.13) 

The entanglement entropy in the deformed theory, (2.11), hence takes a generic form [17]: 

See{A) = + y {{OOH^^) - {OO)) + ... (2.14) 

As we can see, this is only sensitive to the correlators {OOT^u) and {OO). These, however, 
are known to be universal functions in any conformal field theory which only depend on one 
overall normalization constant: 

(0(i,)0fe)) = (2.15) 


and a similar universal function for {0{xi)0{x2)T^y{x^)) whose normalization is also deter¬ 
mined solely by Co [33, 34]. The entanglement entropy of A is hence completely universal up 
to O(A^), only depending on the fundamental constants {a^ , Co} and the operator dimension 
A. 

The divergent structure of the expression (2.14) has been computed in [17] just using 
standard Euclidean GET techniques. However, we must now face an additional subtlety: 
the analysis of [11], which treats the problem at hand using the replica trick and carefully 
analyzes the analytic structure of the resulting partition functions, leads to a different result 
which involves a non-trivial hnite piece. Nevertheless the basic lesson of this more involved 
analysis using the replica trick turns out to be still the same; namely the O(A^) perturbation of 
entanglement entropy is determined solely by correlators {OO) and {OOT^^,). To be precise, 
the hnite contribution to the correction at O(A^) reads [11, 14, 15]:^ 


d+l 


SEEiA) = S^^UA) - A2i?2(rf-A) 


7r^(d-A)r(l + f-A) 

2r(| +h-A) 


+ 


(2.16) 


where R is the radius of the ball A. In this expression the normalization Co has been hxed 
in a way that will prove convenient when comparing with holographic calculations. However, 
we will not need the explicit form of See{A) for the following analysis. 

We will use indices a,b,... for bulk coordinates, .. . for boundary directions, and a, P,. .. for intrinsic 
coordinates on the minimal surface Emin- 

® This is assuming that C? is a relevant operator with A 7^ d/2. See [14, 15] for a discussion of the case 
A = d/2 in holographic theories. 
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2.3 Gravitational explanation for universality 

Having seen that the second order computation of perturbed entanglement entropy of spherical 
regions is only sensitive to {OO) and {OOT^y) correlators (which are universal up to one 
normalization constant), we now want to explain what this means in gravity. 


2.3.1 Holographic parameter matching 

A major hint for understanding the entanglement universality from a holographic point of 
view is the result of [11], where it is shown that for the CFT deformation characterized 
by a coupling A and an operator dimension A as described in §2.2, one can calculate the 
entanglement entropy at O(A^) by solving an auxiliary gravitational problem described by the 
action 


laux — 


^ [fl + - ^ (( 0 ^) 


(2.17) 


The AdSrf+i background geometry gab has a scale which coincides with the cosmological scale: 
^AdS ~ Further, A^ parametrizes the normalization of the scalar action. The prescription is 
to solve first the scalar wave equation in the background gab and then compute its backreaction 
on the metric at next order in perturbation theory to get the change of the area of the 
Ryu-Takayanagi minimal surface. For this task we are therefore not interested in the full 
action (2.17), but only in the linearized equations of motion due to a scalar probe in AdS^+i 
background. In the following all objects with a bar (e.g., Va) refer to the background metric 
9ab’ 

We will now comment on the number of free parameters in this problem and how to deter¬ 
mine them, given that they should induce the CFT parameters at the boundary in a way consis¬ 
tent with AdS/CFT. A priori there are three dimensionless quantities {(^Ads/^p) > ^AdS^^’ 
in the action (2.17). The scalar wave equation in AdS, (□ — = 0, has two asymptotic 

solutions with respective falloff behavior z^+ and z ^~, where the boundary is at z = 0 (in 
Poincare coordinates such as (C.4)) and 


^ ± + 4ds ■ (2-18) 

The boundary condition for the scalar is required to be such that it sources an operator 
with given dimension A and coupling A. By standard AdS/CFT techniques [1], we therefore 
assume® A > d/2, tune the scalar mass such that A+ = A and fix the dominant asymptotics 
as 4>{z —)• 0) ~ A z^~. The mass of the scalar is also immediately fixed by demanding that 
it sources an operator with dimension A; the usual prescription for this is to take —(i^/4 < 
^AdS ^ ^ such that 

el^sm^ = A{A-d). (2.19) 

® The assumption A > d/2 corresponds to performing so-called standard quantization. Note that (d—2)/2 < 
A < d/2 would still be compatible with field theory unitarity. The latter regime can be treated by choosing 
A_ = A. See [13] for a discussion of this case. 
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The remaining two parameters of the theory also need to be fixed by matching with the 
CFT. On the one hand, ^Ads/-^? makes an appearance in the normalization of the universal 
scaling of entanglement in the CFT in terms of the generalization of a-central charge, c.f., (2.5) 
and can hence be fixed by demanding (2.6). On the other hand, parametrizes rescalings 
of the bulk held (p in the same sense in which the dual boundary two-point function (OO) 
depends on an overall normalization Cq as in (2.15). Since p has already been equipped 
with boundary conditions which lead to the correct coupling A for O, we can follow standard 
techniques to work out the normalization in terms of Cq- This involves computing the 
bulk to boundary propagator for p and matching its boundary two-point function with the 
CFT correlator {OO). The result is [1, 13] 


Co 


A - d/2 r(A) 4;^ 
T,d /2 r(A - (i/2) .^^-1 ' 


( 2 . 20 ) 


We stress that, as shown in [11], the auxiliary gravitational problem just described can be 
used for computing perturbed entanglement entropy of ball shaped regions in any CFT. Given 
the universality of the correlators {OO) and {OOT^P) in all CFTs (including those which are 
dual to cosmological Einstein gravity), this result is not very surprising. It is, however, quite 
remarkable that the most natural way of rewriting the CFT problem is in terms of the auxiliary 
gravitational system (2.17). As we shall see, there are many other (higher curvature) models 
of gravity which could equally well serve as the auxiliary system for this problem. Turning this 
logic around, all these different gravitational theories are indistinguishable if only perturbed 
entanglement data of spheres is measured. We now turn to illuminating this point from the 
gravity side. 


2.3.2 Construction of the most general gravitational model 

We now want to make more precise the statement that all physically acceptable gravity models 
are equivalent for the computation at hand. Our approach is to construct a large class of 
auxiliary gravitational systems which make the same predictions as (2.17) in this respect. Let 
us take the CFT data as given and try to construct a holographic description in terms of 
semiclassical gravity ’bottom up’. The first step of such a program could be the question: 
what are the minimal bulk ingredients if they are supposed to compute entanglement entropy 
of spheres in vacuum and after turning on the perturbation? We will now answer this question 
by building general linearized bulk actions which satisfy the following conditions: 

• The field content is a symmetric spin-2 tensor (the metric) and a scalar field, which 
source the boundary operators and O, respectively. The model gives dynamics to 
these helds and defines a coupling between them which allows to compute the linearized 
backreaction hab on an AdS^+i background g^b after perturbing the latter with the scalar 
probe p. 

• There are no ghost degrees of freedom; in particular the linearized equations of mo¬ 
tion on AdSrf+i are second order in derivatives. (See below for a justification of this 
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assumption and §A for generalizations in the case where the higher curvature terms are 
perturbatively small and this assumption can be dropped.) 

• The action is diffeomorphism invariant (c.f., [35, 36]). 

We will find that these conditions specify the linearized action uniquely up to normalization 
constants; see (2.21) for the final result. 

Let us briefly explain why the equations of motion should be second order. It is well known 
that equations of motion which are higher than second order lead to pathological behaviour 
which we want to exclude from our discussion of theories which are dual to unitary CFTs. In 
particular, higher order equations of motion generically give rise to further propagating degrees 
of freedom with negative energy [37-41].^ The presence of ghosts in such theories thus leads to 
classical instabilities [43] (which would lead to negative norm states and therefore non-unitary 
quantum theories upon quantization). These features would also be visible in the dual CFT 
as operators with complex conformal dimension or negative norm. Indeed, via holography the 
problem of perturbed entanglement entropy is described in the CFT by correlation functions 
of the stress tensor T^i, and a relevant scalar operator O. This setup should be described in 
gravity by nothing more than a spin-2 graviton which sources T^i, and a scalar field to source 
O.® If the graviton equations of motion were higher than second order in derivatives, then the 
additional unwanted ghost modes would couple to other (non-unitary) operators. Therefore, 
whatever the full gravitational theory is, after linearization of the action around AdS^+i the 
excitation hab is bound to appear in such a way as to yield second order equations of motion. 

Despite these arguments in favour of two-derivative equations of motion, higher derivative 
models with higher order equations of motion are regularly studied as toy models for semi- 
classical corrections to Einstein gravity induced by some consistent truncation of full unitary 
quantum gravity (which is assumed to be UV complete). Such models can be justified as sen¬ 
sible toy models in a perturbative framework as long as the ghost modes are very heavy and 
do not go on-shell (which turns out to be the case whenever the higher derivative couplings 
are small). We explore this setup in §A, where we show that in this context our arguments 
concerning entanglement entropy still hold. The reason for the latter is that at low energies 
the ghosts can be ignored and the sector of physical modes is exactly the same as for ghost-free 
models which we are now going to analyze. 

^ There are special circumstances where the additional degrees of freedom can be dealt with by using 
appropriate field redefinitions. For instance, f{R) theories of gravity naively have higher order equations of 
motion, but are known to be equivalent under field redefinitions to Einstein gravity coupled to a (positive 
energy) scalar field such that the boundary theory is unitary [42]. The graviton in this case couples not only to 
the boundary stress tensor, but in addition to some scalar operator. The CFT therefore has some more specific 
features than we assumed in the beginning. For validating our conclusions in the most general circumstances, 
it would be interesting to explore the most general requirements for a bulk theory to admit a unitary dual 
CFT. A first step in this direction might be a gauge-fixed analysis where the scalar part of the graviton is set 
to zero, c.f., the transverse traceless gauge employed in §A. This would at least take care of f{R) theories. 

® Of course, it would be possible for O to be some other scalar operator, e.g., one corresponding to the 
square of a scalar field. Since such scenarios would not change our analysis qualitatively, we will assume that 
O is simply sourced by (p with appropriate boundary condition. 
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Let us now construct the minimal ingredients for a linearized action for hab and </>, which 
leads to second order equations of motion. We will argue that the most general action of 
this type is the same as the linearization of Einstein-scalar theory (2.17) up to various nor¬ 
malizations. For illustration, we can hence start by expanding the metric around an AdS^+i 
background solution gab as gab = ffafe + h-afe and linearizing the Einstein action (2.17). Including 
some convenient normalization factors, this procedure leads to® 


I : 


1 

24“^ 






/2 d/2 pd-1 2 ^ 

«AdS ^ ^AdS ^ 

^ /i“^ (-I- [(Vfii)^ -I- rn^d)^] Qah + Vnd'Vhc 


1 ((V0)= 


-|- m 


,2^21 


where Gabcd is the graviton kinetic operator 

Gated 2 i.9a(c9d)b 9ab9cd^ I—^ ^[c9d){a^b) ~ {^9ab'^cSdd “1“ 9cd^b^ 

d (I _ _ - - \ 

“1“ ^£2 ^2 9ab9cd 9a{c9d)b J • 

A few explanations are in order concerning the way we wrote (2.21): 

(i) In writing (2.21) we introduced a normalization \h in the graviton kinetic term for later 
convenience. By expanding the Einstein-scalar action (2.17), one finds simply {\h)Einstein = Ij 
but we will shortly encounter other (higher derivative) theories with other values of \h- 
[a) Similarly, we wrote a constant term proportional to a*^. This is also using the benefit 
of hindsight (in particular (2.7)) to already implement the correct generalization for higher 
derivative theories. For the present case of Einstein gravity, (a"^)Einstein is just given by (2.6). 
For more general theories, notice that the constant term in (2.21) is determined by the require¬ 
ment that switching off hab (p ii^ th® action is equivalent to evaluating the full non-linear 
Lagrangian on the AdS background. By means of (2.7) we know that the constant thus ob¬ 
tained needs to be a)) (with appropriate pre-factors). 

(in) Note also that we did not write the 0{h) term of the action. This term would just define 
the equations of motion for the AdS^+i background geometry gab- By assuming an expansion 
of the form gab = gab + dab, such terms vanish automatically. 


The graviton-scalar action (2.21) was obtained by linearizing the Einstein-scalar system 
(2.17) around a given AdS background with scale ^AdS- However, we will now argue that any 
theory which satisfies the basic conditions listed at the beginning of this subsection, takes this 
form. As we will see, the only freedom is parameterized by the normalizations {a^, Xh, A<^, m^}. 
Of these, a*^ has already been fixed manifestly by using (2.7). Also and -^^dS could be 
fixed in terms of field theory quantities via (2.19) and (2.20). The graviton normalization, 

^ To verify this by hand, the identities in §D are helpful. 
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\h, will be discussed shortly. This will then demonstrate the universality of the dynamics of 
perturbing an AdS geometry with a scalar field up to normalizations. 


Graviton propagator on AdSj^+i: Clearly, the field content in the bulk should involve a 
spin-0 field cj) to probe the geometry, whose boundary condition cl){z —>■ 0) ~ A sources 
the relevant deformation O. But before specifying the matter probe in this much detail, let 
us discuss the general gravitational response to any small perturbation. This is described by 
a spin-2 field propagating on the background gab- We will write this as gab = 9ab + hab, where 
hab is the linearized perturbation. 

We now want to establish the most general physically acceptable equations of motion for 
hab, as described above. It is straightforward to write the most general (up to total derivatives) 
action at 0{h‘^) of the symmetric spin-2 field which involves at most two derivatives. We can 
parametrize it as (see, e.g., [44]) 


I{h^) - 


2e 


T i ^^ahbS^h^^ + ^Vbh'^^Vahl - ^V^hV^hab 




A'P 




- 64 hath' 


ab 


(2.23) 


where h = h°'^gab- All raising and lowering of indices as well as covariant derivatives with a 
bar refer to gab-^^ We can already see that compared to the infinite number of possible higher 
derivative theories of gravity, the above linearized action is quite restrictive, having only six 
free parameters {A/^, 61 , 62 ; ^ 3 ) ^ 4 }- In fact the action will be even more constrained once 

we impose that the full theory should exhibit diffeomorphism invariance. It is a classic result 
that demanding a diffeomorphism invariant theory beyond the linearized level singles out a 
unique graviton propagator [45]. In our explicit notation, we can give a simple justification of 
this statement by considering a spin -2 gauge transformation (aka linearized diffeomorphism) 
hab hab + ‘^d{aCb)- The 0{h?) action transforms nontrivially: 

^hh^) = ii ! {{hi- XhW^Uhab + (61 - h 2 )V aVSdh^<^ + ih - b 2 )Vaah 

P '' (2.24) 

- ( 2 M 2 - Vah + (264 %hi\ , 

V ^AdS / V ^AdS / j 

where we used partial integration and (D.2). In order for the action to be diffeomorphism 
invariant, the above expression has to vanish for arbitrary Since all terms are independent, 
this condition forces upon us five relations among the six parameters: 

hi = b 2 = bs = Xh , = —2 — Xh, 64 = 2 . (2.25) 

^I-Ads 

Note that the coefficient Xh of the first term in (2.23) cannot be vanishing as this would mean that no 
spin-2 degrees of freedom propagate at all. By this we ensure that no gauge choice can eliminate the spin-2 
mode. For instance, we could choose to work in a transverse traceless gauge (Vah“*’ = 0 — h), but the first 
term would still survive as it is the source of the spin-2 component. On a similar note, we need A/i > 0 in 
order for the spin-2 graviton not to be ghost. 


13 








Any other choice of parameters would not correspond to the linearization of a diffeomorphism 
invariant theory. The only free parameter which is left undetermined is the overall scale 
Xh of the graviton action. In a full non-linear theory of gravity this is, of course, fixed by 
the normalization of the action and we will compute some examples in §3 to illustrate this. 
For instance, we have already seen that the graviton action (2.23) is exactly the same as the 
second order term in a perturbative expansion of the Einstein-Hilbert action with cosmological 
constant scale i = .^AdS and overall normalization = 1. 

This construction shows that by just looking at the linearized graviton propagator on 
AdS of any reasonable diffeomorphism invariant theory of gravity, we will never be able to 
distinguish the theory from cosmological Einstein gravity. The underlying reason for this is 
the universality of spin -2 propagation on maximally symmetric backgrounds. 

Coupling to scalar matter: Having constructed the most general quadratic action for the 
metric perturbation, we should now in a similar fashion argue that ( 2 . 21 ) couples hab to the 
scalar field in the most general way at this order. There are several possible Yukawa couplings 
and self-interaction terms. However, for consistency of our perturbative scheme, the only 
terms we are interested in at this order are (j)(j) self-interactions and /i(/>(/> couplings. We will 
now describe the most general actions of these kinds: 

• Scalar field kinetic and mass term: At lowest order in perturbation theory (the 
probe limit), we are only concerned with couplings of cfxf) to the background geometry. 
Up to total derivatives, the most general such action with at most two time derivatives 
may be parametrized in terms of a mass parameter m? and an overall dimensionless 
normalization A^: 


V) = ^ [-1 } • ( 2 - 26 ) 

This action defines the background dynamics of </> as (□ — m?)(j) = 0. 

• Backreaction on the geometry: There are several possible couplings of the type 
hijxl), which are inequivalent upon integration by parts. We can parametrize them as 
follows: 

= Ail [ I - 7(ci(V(())^ C2Tr{^4>‘^) -h y h'^^Va^Vb^ 

^ (2.27) 

-I- C4 h°-^(j)Va'^b4> + C5 , 

Clearly the coefficients C 2 and C 5 are on-shell equivalent. We will thus set C 5 = 0 
in what follows. Now demanding diffeomorphism invariance, an argument completely 
analogous to the one we gave for (2.23) implies the following conditions for the remaining 
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parameters: 


Cl — C2 — C3 , C4 — 0 . 


(2.28) 


which leaves us with only one free parameter ci. However, it is immediately clear that 
the overall normalization has to be the same as in (2.26), i.e.. 


Cl = . 


(2.29) 


If this normalization was any different, hab would not be the response to a scalar with 
mass m? and normalization A^. Taking into account these constraints, we can rewrite 
(2.27) as 


with T<J> 


^ f I —H f,ab 

[(V()))^ + 5a6 +Va(/>Vfe(/>, 


(2.30) 


Summarizing these two points, we recognize the most general scalar-graviton couplings 
which are consistent with the perturbative scheme and diffeomorphism invariance as being 
exactly the same as the expansion of the following action up to 0{h): 


I matter — + 0(h ) 


1 




-i((V0)2 + mV) 


+ 0(/i2 


(2.31) 


where gab = gab + hab- Of course, this result is not too surprising as this is the standard action 
for a massive scalar minimally coupled to gravity. 


Most general second order action: We can summarize the above discussion as follows. 
The most general way of coupling hab to a scalar (p in a way that is consistent with our pertur¬ 
bative scheme and respects diffeomorphism invariance and some very basic constraints coming 
from the particular CFT computation we are interested in, is given by (2.21). This action 
has only four dimensionless parameters {^Ads/'^P j ^^-^AdS ^ ^(t> j^h}- We already explained in 
§2.3.1 how for our problem the first three of these have to be determined holographically in 
terms of the dimensionless CFT quantities {a^, A , Co}- We will argue shortly that the fourth 
parameter \h does not enter the computation of interest and therefore its value is irrelevant 
for computing perturbed entanglement entropy. In addition, there is, of course, the boundary 
coupling A which enters the gravitational setup by prescribing the boundary condition for (p. 

As a consequence, any covariant theory of gravity whose linearization takes the form 
(2.21), provides a good auxiliary system for computing SEEi-^) at second order in perturbation 
theory. For instance, the action (2.21) with A/j = 1 is precisely the same as the linearization 
of the cosmological Einstein-scalar action (2.17) via gab = (jab + hab- There are, of course, 

Alternatively, one can determine the coefficients by demanding conservation of the matter stress tensor, 
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many other non-linear actions whose linearization on maximally symmetric backgrounds also 
takes the form (2.21) with a different value of e.g., those of Lovelock type [46] or quasi- 
topological gravity [47]. Indeed, even if the CFT in question is actually dual to some particular 
higher derivative theory; for the purpose of calculating the perturbed entanglement entropy, 
one can always work with the Einstein action (2.17) after suitably matching the parameters 
{^Ads/^Pj ^-^^AdS ’ with the CFT boundary conditions. The different linearized actions 
will then only differ by their value of Xh- We will see examples of this in §3. 

Since the normalization X^ is the only ingredient in (2.21) which would genuinely dis¬ 
tinguish bulk theories with higher derivative couplings from those without, we now need to 
understand why the linearized computation of perturbed entanglement entropy of a ball is 
not sensitive to this. 


Holographic interpretation of A^: Before we continue with the gravitational analysis, it 
is instructive to elaborate on the meaning of X^, which has a very natural interpretation from 
the CFT point of view. The graviton normalization X^ corresponds to the normalization of 
the CFT stress tensor two-point function, usually denoted by Cj,. It sets the strength of the 
leading divergence in the two-point function of CFT stress tensors: 

{T^^{x)Tp^{0)) = ^Ipu,pa{x), (2.32) 

with some well known universal function Tpy^p„{x) [33, 34]. In d = 2, the normalization 
agrees with the usual c-central charge, = c. Similarly, in d = 4 one finds = (40/7r^)c. 
Once a holographic description of the CFT has been given, usual AdS/CFT methods can be 
used to compute the left hand side of (2.32) as a bulk graviton correlator. Since Xh sets the 
normalization of the graviton propagator in the presence of matter sources (c.f., (2.38)), it is 
rather clear that C^, will be proportional to Xh- Doing the holographic calculation in detail 
leads to the standard result [48] 


= 


2 e 


nd-l 

'^AdS 

id—1 


Xh , where 


^ 2d(d + l) F(d) 

“ d - 1 7r'^/2r(d/2) 


(2.33) 


It is now obvious that once ^Ads/(■'P has been fixed such as to reproduce the other central 
charge via matching (2.7), there is no free parameter left in (2.33) which could absorb the 
normalization Xh- In our setup, gravity theories with different Xh will hence describe CFTs 
with different Cj,- Conversely, a CFT measurement of Cj, would allow to for a distinction 
between different bulk models. However, we will now argue that the perturbed entanglement 
entropy of spheres is agnostic about this. 


2.3.3 Second order Wald entropy is independent of graviton normalization 

As we have seen, any gravity theory whose action linearized around AdS^+i takes the form 
(2.21) can be used to compute second order perturbations of entanglement entropy of spheres. 
However, we identified one parameter Xh in the action (2.21) which sets the graviton normal¬ 
ization and turns out to be sensitive to more detailed bulk structures such as higher curvature 
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terms. We will now argue that this parameter does not enter the second order entanglement 
entropy.^^ This is the gravity version of the statement that (2.16) is independent of C^,. 

We recall that due to the spherical symmetry of A, extrinsic curvatures on the minimal 
surface £min vanish and the entanglement entropy functional to be evaluated on £min is just the 
Wald functional. In vacuum this reduces to just the area of £min multiplied by a normalization 
which is basically determined by as in (2.5). The perturbation can then be computed in 
holography as 

^SEEiA) — 5 Swaldi£min) — ^ i ^TT f \/T Ttj ^ab^cd\ • (2.34) 

Clearly 5Swald{£min) can be expanded perturbatively by considering different powers of hab 
and (p. From the structure of the gravitational backreaction it follows that there is no con¬ 
tribution at 0(A). At O(A^) (i.e., at first order in hat) the variation in (2.34) acts on the 
objects inside the integral, but it does not change the surface £min itself.^^ Since the surface 
£min is maximally symmetric, we anticipate that this variation of the Wald functional is just 
a variation of the area form (up to normalization). That is, the evaluation of (2.34) should 
yield precisely the same answer as in Einstein gravity up to normalization. Let us make this 
slightly more explicit to figure out the correct normalization. 

To evaluate (2.34) at O(A^), we should expand the integrand linearly in the perturbation: 


SSwald{£min) — 27r 




6C 

—— 7T,(jfe IT'cd 

0£iabcd 


0(h) 


(2.35) 


The obvious way to evaluate this is to start with the full non-linear Lagrangian, expand the 
Wald entropy density around the AdS background, and then identify the 0{h) term in (2.35). 
However, even without knowledge of the full theory it is clear that ultimately the integrand 
can only be sensitive to the part of C which is at most quadratic in hab, he., we should be 
able to predict the answer just based on the most general second order Lagrangian (2.21). 
Moreover, only the purely gravitational sector of (2.21) is relevant for Wald entropy at O(A^). 
This is quite clear from the observations that a minimally coupled scalar field in gravity does 
not couple to curvature and the action (2.21) agrees with that of a minimally coupled scalar 
at the order we are interested in. 

From these observations, we are led to conclude that the integrand in (2.35) formally has 
the same structure in any allowed theory of gravity up to the normalization Therefore, 
without loss of generality, we can compute the integral as if it came from Einstein’s theory 
and multiply by A^. But in that case, the answer for 5Swald{£min) is well known: it is just the 
variation of the area functional. Therefore, (2.35) can be rewritten in general as a variation 

See §C for related observations based on the entanglement first law. 

Variations of the bifurcation surface only affect the Wald entropy at the next order, which is beyond our 
discussion. To see this, note that the surface Emin extremizes the Wald entropy functional upon variation with 
respect to the metric. First order variations of the metric (which are O(A^) in our conventions) hence leave 
Emin unchanged. 
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of area, multiplied by the normalization Xh that encodes the higher derivative dependence: 


SSwaldiSrmn) = 27r • (2-36) 

V / ^AdS 

where 7 q,^ is the background metric induced on Emin by gab and is the induced metric 
perturbation^^ Note that the arguments just given were heuristic and we are discarding 
boundary terms which would be important to regulate the divergences. A detailed derivation 
of (2.36) can be found in §B. 

The integral (2.36) can now be used to derive the independence of SSwaidiEmin) on Xh- 
We just need to find the on-shell dependence of h^p on X^ and combine it with the prefactor in 
(2.36). To this end, consider the linearized equations of motion implied by the general action 
( 2 . 21 ): 

\hQab^^hcd = -^T^'’ . (2-38) 

From this equation it is clear how to go from the universality of graviton propagation to 
the universality of perturbed Wald entropy: apart from the overall Xh every dimensionless 
parameter in (2.38) has been fixed by the boundary input {a^,A,Cc)}. The combination 
Xhhcd therefore satisfies a universal equation of motion which is independent of the non-linear 
details of the bulk theory. Since the perturbed Wald entropy (2.36) depends on precisely 
this particular combination it is clear that SS\Yald{Emin) will be the same in all theories we 
consider. Said differently, the function hab which satisfies the backreaction equation (2.38) 
is proportional to 1/Xh such that the Xh dependence in (2.36) cancels. We conclude that 
5See{-A) = ^Sivald(Emin) is independent of Xh- This is consistent with the field theory result 
(2.16) being independent of C^,. 

This concludes our analysis of the role which various dimensionless parameters play. Hav¬ 
ing verified the cancellation of Xh, the explicit solution of (2.38) and evaluation of (2.36) can 
be done without loss of generality at = 1, i.e., Einstein gravity. We refrain from repeating 
this calculation. See, for example, [14, 15] for an explicit evaluation of (2.36) in the case of a 
scalar field source. The result of their computation perfectly matches the CFT answer (2.16). 


3 Example: Gauss-Bonnet theory 


In this section we wish to illustrate the general argument, using the example of gravity theories 
with Lagrangians quadratic in curvature. This discussion will make contact with a number of 
previous higher derivative computations of entanglement entropy. In particular we will show 

Explicitly, if denotes the coordinates of Emin embedded in AdSd+i, we have 


_ 


For more conventions see §B. 


_ _ 

9ab , 


h, 


(2.37) 





that various results for perturbations of entanglement entropy of spherical regions dual to 
higher derivative theories can be equally well reproduced from a calculation in cosmological 
Einstein gravity by renormalizing the discrete set of free parameters which play a role at second 
order in perturbation theory, i.e., {^Ads/^P j ^-^-^AdS ’ appearance of 

the additional parameter Xh which drops out of the formula for entropy. 


3.1 Two gravity models 


The basic gravitational model that we want to take as the reference system is just the auxiliary 
Einstein-Hilbert theory described by 


I = 


1 

2i^-^ 




R + 


d{d - 1) 


Y (( V(/>)2 0 ^) 


(3.1) 


with asymptotic boundary condition (j){z —)• 0) ~ A z'^~^ such that (p sources a scalar operator 
with dimension A and coupling A in the GET. We parametrize the lowest order backreaction of 
the scalar probe on the geometry as gab = Qab + hab, where (jab is AdS^+i with scale ^^dS ~ 
The independent dimensionless parameters in the gravitational problem can therefore be taken 
as {(-^Ads/^p) 5 "^^^AdS ’ while A/i = 1 for Einstein gravity. 

In this section we wish to reproduce predictions of four-derivative theories by using the 
above Einstein theory with renormalized parameters. The only four-derivative theory which 
passes the consistency condition of yielding second order equations of motion is Gauss-Bonnet 
theory, which we will discuss here. Nevertheless, in §A.l we explore the most general four- 
derivative theory of gravity and come to very similar conclusions in the case where the higher 
derivative terms are small perturbative corrections to (3.1). 

With this motivation in mind, let us take the following action as the theory to compare 
with: 


' 24-^ / ^ 

where C^.d. = 7 - ^RabR'^^ + RabcdR''^'''^) 


R + Cb.d. - ^ {{Vpf + (t>^) 


(3.2) 


again with a fixed boundary condition (j){z ^ 0) ~ Xz‘^~^. The AdS^+i background solution 
of this theory will have some other scale ^AdS which is a function of i and 7. In the context of 
both theories I and I, we will use gab to refer to the AdS^+i background. From the context 
it should be clear whether (jab is the background with .^AdS or with ^AdS- 

Our basic aim is to show that for the purpose of carrying out the perturbed entanglement 
entropy computation of spherical regions, the actions (3.1) and (3.2) are indistinguishable if 
the dimensionless parameters {(^Ads/^p) ■, rrr^^dS ’ suitably expressed as functions of 

{(^Ads/^p)) , A^ , 7} such that they are consistent with the GET boundary conditions. 

In order to show this, we will now derive the relevant ingredients for such a computation order 
by order in A from both I and I. We will demonstrate that the parameters of Einstein gravity 
can always be chosen such that the two perturbative expansions match. The higher derivative 
coupling 7 can be absorbed entirely into renormalizations of the other dimensionless quantities. 
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3.2 Parameter matching 

We will now determine the relation between the dimensionless parameters of the models I and 
I such that both make identical predictions for vacuum and perturbed entanglement entropy 
of spheres. 


Finding the AdS^+i backgronnd solution: Before turning on the scalar perturbation, we 
need to determine the maximally symmetric AdS^+i background solutions of the two theories 
(3.1) and (3.2). It is well known that such solutions exist in both theories and are ghost-free 
[49] if we choose the respective AdS scales 



1 +v'l-4(d-2)(d-3)7) . 


(3.3) 


We will henceforth expand around the AdS^+i background and thus treat .^AdS i as 
equivalent (and similarly with tildes). 


Matching of Wald entropy (or a[)): Let us remind the reader of the result of §1 that Wald 
entropy of maximally symmetric spaces is always proportional to area. We can thus compute 
the unperturbed entanglement entropy of A as the Wald entropy of hyperbolic black holes 
corresponding to the two AdS radii (3.3) of Einstein and Gauss-Bonnet gravity, respectively. 
By integrating (2.5) for the two gravity models, one finds 


\ 

2(d-l)(d-2)7^ . 

^AdS/ 

(3.4) 


Q * Vol(]Hd_i) 

^Wald = —(“/4) 


27r 


nd-1 

"^AdS 


Swald = ^r{d/ 2 ) d*. 


TT' 


Vol(Hrf-i 

nd-l 

^AdS 


— 


— 


TT 


d/2 


pd-1 

^AdS 


r(d/2) 4-1 ’ 


TT 


d/2 


pd-1 

^AdS 


r(d/2) 


1 - 


Given the entanglement entropy S'gg(A) with fixed normalization in the GET, it can be 
reproduced by the computation of Wald entropy in either of the gravity models by requiring 


pd-l 

^AdS 

pd-l 


pd-l 

•^AdS 


pd-l 


^l-2(d-l)(d-2)7 



(3.5) 


This determines the dimensionless ratio (-Ads/^'P in Einstein gravity in terms of {^Ads/^P j 7} 
in Gauss-Bonnet theory such that both theories make the same predictions for Wald entropy 
of the maximally symmetric horizon. From the GET point of view, this matching is the 
condition that is a fixed fundamental parameter of the GET which should be reproduced 
by either candidate of a gravity dual. 
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Matching of scalar coupling (or Cq)- Perturbing the CFT with O means introducing 
a massive scalar in the gravitational action. The mass of the scalar is constrained by the 
requirement that the boundary value of the scalar sources the operator O with given dimension 
A. This fixes 

4ds -d)= 4ds ■ (3-6) 

The dynamics of the scalar field is derived at lowest order from the propagation on AdS(i_|_i 
vacuum: 

(□ — m^) (/> = 0 , ^ 0) ~ A 2 '^“^ , (3.7) 

and similarly with tildes. Once the boundary condition is imposed in this way, standard 
holographic techniques lead to the matching between the normalization Cq of the relevant 
CFT correlators and the normalization of the scalar field action via (2.20): the two gravity 
theories I and I lead to the same Cq provided that 


pd-l pd-l 

\ ^AdS _ \ ^AdS 

pd-l pd-l 

-v-D •v-'D 


(3.8) 


so the ratio should be fixed to be the same as in (3.5). 


Computation of perturbed entanglement entropy: We finally want to check explicitly 
the arguments of §2.3.3 and verify that after the matching conditions above have been imple¬ 
mented, second order entanglement entropy is automatically the same for I and I. We will 
also verify that the difference between A/i = 1 and \h cannot be absorbed in any further free 
parameters. This will lead to the conclusion that the central charge Cj, is genuinely different 
in the CFTs dual to I compared to those dual to I after the above matchings have been 
performed. 

Our perturbative scheme in the CFT demands that we consider the linearized equations 
of motion in the bulk. The linearized equations of motion can be obtained from the actions 
(3.1), (3.2) by linearizing around We find for I and I, respectively [50]: 

XhGab"Xd = -^T^f , 

( (3.9) 

\hQab^‘^h,d = —^ f^\ 

where and are the source terms (2.30) for scalar matter with mass and fh?, 
respectively propagating on AdS^+i with scale ^AdS ^md ^AdS- The kinetic operators Qabcd and 
Gabcd are defined as in (2.22), but, of course, also refer to backgrounds with respective scales 
^AdS and ^AdS- The coefficients \h and \h in (3.9) can be computed for I and I theories and 
turn out to be 

Xh = l, Xh=(l-2{d-2){d-i)^^\ . (3.10) 

V 4ds/ 
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Clearly there is no free parameter left in (3.9) which could be renormalized to absorb these 
normalizations. From (2.33) it is hence clear that the C^, charges of the dual CFTs will 
differ by the ratio X^/Xh- At this stage, any CFT observable which did depend on Cj, would 
therefore genuinely distinguish between the two gravity models. As shown before, 6See{-^) 
does not depend on Cj,, though. 

Our matching prescriptions (3.3), (3.5) and the graviton kinetic term scaling (3.10) have, 
of course, appeared in the literature before. For example, see [20] for a computation of 
perturbed sphere entanglement entropy in Gauss-Bonnet theory. One can easily verify that 
their computation gives an answer which takes the form of perturbed entanglement entropy in 
Einstein gravity with dimensionless parameters rescaled as predicted by our general analysis.^® 
At this order the entanglement entropy computation for spherical regions is therefore not 
sensitive to the Gauss-Bonnet coupling. 

4 Discussion 

Let us summarize our arguments and comment on implications for the program of bulk re¬ 
construction. 

Summary: We considered the second order perturbation of GFT vacuum entanglement 
entropy of spheres due to scalar operator deformations. This quantity exhibits a universal 
functional form which only depends on a small number of fundamental GFT parameters to 
set various normalizations. We explored a similar universality statement about gravity, which 
explains why the corresponding holographic computation can always be recast in terms of one 
universal gravitational system, viz., linearized Einstein gravity coupled to a free massive scalar 
field. We argued that any ghost-free gravitational theory other than Einstein gravity can be 
used to make identical predictions due to the uniqueness of diffeomorphism invariant spin-2 
propagation on AdS. For instance, we demonstrated how in Gauss-Bonnet theory both the 
vacuum and the perturbed sphere entanglement entropy (obtained by computing Wald entropy 
of a minimal surface) can be reproduced from just Einstein gravity by suitably renormalizing 
the dimensionless parameters of the theory. In this sense, knowledge of this entanglement 
data does not allow to distinguish between Gauss-Bonnet theory and Einstein gravity. 

More precisely, the perturbed sphere entanglement only depends on the dimension A of 
the deforming operator O, the generalized o-central charge (in vacuum) and the normal¬ 
ization Co of the (OO) correlator (at second order in the perturbation). We saw that any 
viable gravity model has enough dimensionless parameters to account for these dependen¬ 
cies (the scalar mass the normalization of the gravitational action, iAds/^v-, and the 

normalization of the scalar field, A^, respectively). In particular, no higher curvature terms 
in the bulk are required in order to obtain a sufficient number of couplings for consistently 
performing the holographic matching. In this problem, the only signature of higher curvature 

See, e.g., Eq. (7.7) of [20] with the understanding that their Gauss-Bonnet coupling A equals 2'y in our 
notation. 
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couplings would be a renormalization of the parameters of cosmological Einstein gravity. This 
implies that perturbed entanglement entropy of spheres can be computed using any ghost-free 
gravitational model as an auxiliary system, even for CFTs without any simple gravity dual at 
all [11], 

In face of this, it is clear that for the derivation of bulk dynamics from entanglement in the 
sense of AdS /GET it is desirable to reconstruct bulk features which are not entirely universal 
in the above sense. The obvious question arises which quantities one should consider instead 
of or in addition to entanglement entropy of spheres in order to achieve a dependence on bulk 
features which are not the same for all CFTs. We will come back to this question below. 

Generalizations: There are some modifications of our analysis which would be interesting 
to explore. For instance, one could compute the perturbation to higher order in A. By 
extending the perturbation theory of §2.2, it is clear that entanglement entropy will then 
be sensitive to higher point correlation functions which would probe bulk dynamics beyond 
the linear approximation. While this would be very interesting to see explicitly, the major 
complication in the bulk would be a backreaction not just on the Wald entropy density, but on 
the minimal surface itself. Once the minimal surface has non-vanishing extrinsic curvatures, 
ambiguities in the entanglement entropy functional need to be addressed [21-24]. By resolving 
this issue, one may hope to derive bulk dynamics beyond linear order. 

We mainly focused on perturbations of entanglement due to deformations of the CFT by a 
relevant scalar operator. As another generalization, one can imagine deformations by operators 
of other spin. As long as the deformation does not involve a stress tensor, our conclusions are 
expected to hold since the new operator comes at least with an overall normalization which 
is all that is required to account for the renormalization of .^Ads/-^P- Consider, however, the 
case where the CFT is deformed by a stress tensor insertion. In that case, the first interesting 
perturbation of the entanglement entropy occurs at 0(A) instead of O(A^), where one gets a 
contribution ~ (TT) with normalization set by Cj,. In the bulk the problem of gravitational 
backreaction changes in a similar way and the extremal surface will feel a perturbation at 
0(A). It would be interesting to explore this case in detail and consider the dependence on Cj, 
to draw similar conclusions about the effect of higher curvature couplings. The dependence on 
Cj, is expected to lead to a way of distinguishing Einstein gravity duals from higher curvature 
ones. In fact it has been argued before (as we review in §C) that generic perturbations of 
spherical region entanglement entropy can be used to derive linearized bulk equations including 
their normalization set by Cj [6]. 

Further generalizations could involve perturbations of the CFT other than operator de¬ 
formations. For instance, a small geometric perturbation of the spatial sphere should lead 
to similar conclusions as those just mentioned. Similarly, 1/N corrections have recently been 
explored in order to improve on the entanglement first law [51]. It would be interesting to 
pursue this further and explore ways to reconstruct bulk dynamics at a level where one can 
distinguish between different higher curvature theories. Instead of considering individual ap¬ 
proaches in detail, let us now outline in general how a good choice of CFT observables can be 
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made which are sensitive to interesting bulk features in a controlled way. 


How to “bootstrap” bulk theories: As we have just seen, there are, of course, simple 
ways to circumvent the problem of perturbed sphere entanglement being insensitive to the 
desired detailed bulk information. Roughly speaking, one needs to consider more detailed 
CFT structures in order to distinguish between different bulk models. In general one could 
imagine bootstrapping the bulk theory using the CFT input. Clearly it is hard to derive bulk 
dynamics directly from fundamental CFT parameters such as two-point functions and OPE 
coefficients. But entanglement entropy seems to provide a physically meaningful quantity 
with well-understood universal properties on both sides of the duality, which depends on 
these fundamental building blocks of the CFT in a useful way. At lowest order, the vacuum 
entanglement entropy of spheres can be used to measure the generalized a-charge a^. In our 
particular context we used this to fix the normalization ^Ads/^P of any possible bulk dual 
(whose dynamics is uniquely fixed in the same way as the functional form of entanglement 
entropy is universal). As we have seen, a small perturbation of this entanglement entropy 
does not affect the dynamics enough to determine if there are any higher curvature terms 
present in the bulk theory. To answer this question, one would have to measure some other 
entanglement entropy which is sensitive to further details of the bulk dynamics. 

One natural CFT parameter that we have only encountered briefly so far is the stress 
tensor central charge. Let us elaborate how any quantity which is sensitive to Cj, will also be 
sensitive to bulk higher derivative couplings and hence will not be reproducible by just Einstein 
gravity once ^Ads/-^? has been fixed by matching vacuum entanglement entropy of spheres as 
outlined in this note. In our context, the ratio computed holographically is independent 

of -^Ads/^P cind genuinely takes different values in the presence of higher curvature terms. In 
fact it is clear from (2.7) and (2.33) that higher curvature terms in the bulk are necessary 
whenever jVidjl^fd does not hold. This is, of course, just a generalization 

of the well-known statement that higher curvature terms in the bulk theory are required 
whenever the a- and c-central charges parametrizing the boundary conformal anomaly differ 
[52-54]. 

Any set of CFT observables which depends on both and Cj, can hence be used to 
detect if higher curvature terms in the bulk are necessary. The problem of genuine bulk 
reconstruction is therefore not at all the restriction to quantities which are universal in the 
sense that they only depend on an overall normalization. But it is now clear that one has to 
consider enough such quantities in order to see a necessity to account for their values in terms 
of higher curvature interactions. 

There are many examples of entanglement-related CFT quantities which one may envisage 
to compute in this context. For instance, the entanglement entropies of other surfaces such as 
infinite strips or cylinders are sensitive to both and Cj, [19, 20, 55]. Further, the presence 
of corners and kinks in the spatial region of interest leads to divergences which exhibit a 
universal dependence on the opening angle and (at least in holographic CFTs) depend on 
and Cj, in an interesting way [56-63]. Other interesting quantities with universal features and 
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dependence on Cj. are Renyi entropies of spherical regions (e.g., in the entanglement entropy 
limit g —)• 1) [64-66]. All these quantities have in common a contribution with a universal 
functional form and a dependence on few CFT parameters (in particular and C^). They 
could hence be used to distinguish between CFTs with a dual description in terms of Einstein 
gravity and CFTs which require higher curvature terms in the bulk. 

However, note that even if some CFT measurement of and Cj, leads to the conclusion 
that there can exist a semiclassical bulk theory only if contains higher curvature interactions 
are present, it is usually not possible to determine precisely what higher curvature terms are 
required. Indeed, the signature of, say, a Gauss-Bonnet term on can equally well be 

reproduced by some other higher derivative term. For illustration of this statement see also §B, 
where it is shown how any higher curvature theory simplifies to the background field expansion 
of Gauss-Bonnet form if one is only concerned with linearized equations of motion. In this 
sense Gauss-Bonnet theory is the most general theory which can account for arbitrary values 
of and Cj, in the CFT. In order to further bootstrap the precise bulk dynamics, one needs 
again access to more CFT data (e.g., entanglement data which depends on further 2-point 
functions or OPE coefficients). The reason is, again, that gravitons propagating on AdS have 
very universal dynamics. This is a manifestation of the familiar problem that entanglement 
entropies which are computable tend to be very universal and contain only discrete information 
about the actual bulk dynamics. 

It would be interesting to use this philosophy to derive more detailed (non-linear) bulk 
dynamics for some simple CFT with known holographic dual from first principles. Generalizing 
this approach may then ultimately lead to a genuine reconstruction of the bulk geometry and 
its dynamics by using a suitable set of entanglement entropies. 
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A Generalization to bulk theories with higher order equations of motion 

In the main text, we considered the case where the equations of motion for the metric per¬ 
turbation in the bulk were strictly second order. The motivation for doing so was to exclude 
ghost modes which would occur if the graviton equations were higher than second order. How¬ 
ever, if the higher curvature contributions to gravity are parametrically small, they can be 
interpreted, e.g., as toy models for perturbative string theory corrections to Einstein gravity. 
The ghost modes are then an artefact of the low energy truncation and can be ignored, which 
allows for a considerable increases of the number of acceptable gravitational models covered 
by our analysis. In order to illustrate this, we first generalize in §A.l the Gauss-Bonnet theory 
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of §3 to the case of all Lagrangians which are quadratic in curvature. In §A.2 we then give a 
more abstract argument to show that this example actually illustrates the general case. 


A.l Example: general curvature squared theories 

In this section we generalize the analysis of §3 to theories where the gravitational sector is 
described by 


where Ch.d. = aR^ + ^ RabR^^ + 1 {R^ - + RabcdR^'^" 


(A.l) 


We do this to explicitly demonstrate that such theories - despite the presence of massive ghost 
modes - are acceptable in a perturbative sense. Indeed, we will see how the ghost degrees of 
freedom are suppressed and hence can be neglected for the computations we want to perform. 

We want to model the theory (A.l) with the cosmological Einstein-scalar theory (3.1). 
That is, we will give the matching conditions which express the parameters of the latter 
theory, {(^Ads/-^p),“ terms of {(^Ads/^p), hl^^AdS > ~^h , A</„a,/3,7} such that 
the theories make the same predictions for entanglement entropy of spheres in vacuum even 
after deforming the theory with a scalar operator. The main new ingredient is the presence 
of four-derivative terms in the equations of motion. 

Before discussing the linearized equations of motion, let us give the analogue of the match¬ 
ing conditions (3.3), (3.5). The stable AdSj^+i background solutions in the two theories have 
the following scales in terms of the cosmological constant: 


^2 _ p2 

^AdS — ^ 


^Ads —(^TT^y (d(d + 1 ) a + d /3 -b (d - 1 ) (d - 2 ) 7 ) 


(A.2) 


The generalization of the matching condition (3.5) then becomes 


^ = ^ (l-2(d(d + l)a + d/3 + (d-l)(d-2)7) 


p2 

^AdS, 


(A.3) 


Exactly as in §3.2, the previous equation also immediately fixes the ratio of scalar couplings, 
i.e., is given by the value bracket on the right hand side of (A.3). This guarantees that 

holographic computations of {OO) yield a unique normalization Cq- 

Consider now the linearized graviton equations in the two theories. For this, we hnd it 
convenient to gauge fix the metric. In the transverse traceless gauge (Va/i“^ = 0 = h) the 
universal linearized equations of motion (3.9) now generalize to the following two equations 
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for Einstein and curvature squared theories, respectively [50]: 


Ah □ + 


^AdS 


Ah □ + 


^AdS. 


hab + P n + 


/2 

^AdS, 




ha, = - A* fW . 


(A4) 


where the coefficients Ah and Ah in (A.4) are now given as 

Ah = 1, Ah = I 1 — + l)a + (i/3 + (d — 2)((i — 3) 7 )—j . (A.5) 

We now want to argue that the second term in the second line of (A.4) can be neglected for 
our computation of entanglement entropy. Indeed, this term is obviously suppressed if (5 is 
treated as a small parameter and in the physical on-shell limit (□ -|- )/iah —>• 0 (see also 

[61]). Further, we can also rewrite the second line of (A.4) in the following factorized form: 


From this expression one can anticipate the two poles of the graviton propagator. The Fierz- 
Pauli type mass term ~ corresponds to the physical graviton mode. The other mass 

term ~ (2/^AdS ^h /schematically leads to the ghost which is very heavy for small /3 

and hence does not go on-shell at low energies, ft can be seen as an artefact of truncating 
some well-behaved theory at four-derivative order. 


A.2 General argument 

Despite the previous subsection appearing to be restricted to the specific case of curvature 
squared theories, the conclusion is actually more general, ff higher curvature terms, which 
appear in the gravitational action with some dimensionless coupling a, lead to ghosts modes, 
then the latter generically have a mass which scales with 1 /a. If these modes are interpreted 
as due to consistently truncating perturbative quantum gravity (i.e., |a| <C 1 ), the ghosts are 
unphysical and indeed never go on-shell at low energies. In such a perturbative framework the 
presence of ghosts is hence not due to an instability of the theory (the underlying full theory 
is assumed to be unitary), but it is just a rather harmless consequence of neglecting certain 
sectors of the theory which are not being probed at low energies. 

To illustrate these statements, consider, for example, a gravity sector of the form Cgrav = 
^gravidab, Rabcd) such that the equations of motion are fourth order in derivatives. The most 
general four-derivative gravitational equations of motion linearized around AdS are precisely 
the same as those captured by the curvature squared theories in (A.l). This can be shown 
using the method of background field expansion; see, e.g., [67], where it was indeed argued 
that the linearized action of any higher derivative theory of this type is the same as that of 
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the most general curvature squared theory.^® Concretely, given some such higher curvature 
theory, the action expanded to second order around AdS takes the form (A.l) with parameters 
{a,/3,7} determined as linear combinations of the original higher derivative couplings. We 
can therefore ignore the ghosts of these theories as long as the higher derivative couplings are 
viewed as a semiclassical truncation stemming from some unitary UV-complete theory. 


B Perturbations of Wald entropy: background field expansion 


This appendix complements the analysis of §2.3.3 by computing the expansion of the Wald 
entropy density at 0{hab), derived from the general second order Lagrangian (2.21). This 
analysis will also serve to illuminate how the universality of graviton propagation on AdS 
simplihes the method of background field expansion, thus demonstrating a more general uni¬ 
versality statement about quantities which depend on both and Indeed, we will show 
in what sense Einstein-Gauss-Bonnet theory is the most general theory at Oih?). Any other 
physically acceptable higher curvature theory has the same background held expansion at this 
order. 

It is straightforward that we can focus on the pure gravity part of the general Lagrangian 
(2.21) since the minimally coupled matter sector does not couple to curvature. We wish 
to give a similar argument for the constant term. Ignoring the matter sector puts us in 
the realm of [68] (see also [6]), where a detailed analysis of perturbed Wald entropy on a 
maximally symmetric AdS background has been performed. The best strategy to compute 
the perturbative expansion of Wald entropy is to rewrite the expansion in hat as a background 
held expansion of curvatures Rated around —j^igacObd ~ 9ad9bc)- More explicitly, instead of 

^AdS 

writing explicit metric perturbations hat we write the Lagrangian solely in terms of 


^Rabcd — Rabcd 


f2 

^AdS 


{,9ac9bd 9ad9bc) 


= 


abed T ^2 {dac hbd T 9bdhae 9adhbe 9behad) T ) , 
^AdS 


(B.l) 


where 9ab = gab + hab is the full metric and indices of expressions like this are raised with the 
full metric. This way, the pure gravity part of our general quadratic Lagrangian (2.21) can be 
interpreted as an expansion in ARabcd and its contractions and derivatives. To wit, we can 
always write the pure gravitational couplings of (2.21) as a linear combination of all possible 
combinations of ARabed which contain terms relevant at 0(/i^). Clearly such combinations 
contain at most two factors of ARabed, but they may contain any number of derivatives in 


In this context, see also §B. 
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general: 


V™ - 2^J_, ] 'Ts “ 2 '‘c<i + 0(ft )j 

= ^/^(4']i+4'2)' 

•^ 1 , f ^AiJ = 00 + 02 Ai? 
with < . 

= bi 4ds + b2 ilds + b3 ilds AR^’^^^ARabcd + ... + Oih^) 

(B.2) 

Let us explain this way of parameterizing the action piece by piece. First observe that from 
the form of (B.2) it follows that Oq is just 


oo — — 


2d r(d/2)4~^ 

4ds <dS 


Od' 


(B.3) 


Similarly, one can easily verify that the vanishing of the linearized Lagrangian (i.e. the equation 
of motion for the background solution gab) demands oq = —{2d/£\^g)a2- This fixes the form 
of and shows that the latter is all that would ever show up in the case of Einstein 

gravity (with an overall pre-factor 02 ). Therefore, is the term which parameterizes 

contributions to the graviton propagator which only come from higher curvature terms. More 
tensor structures involving derivatives and contractions of ARabcdARefgh could be written in 
(B.2) in order to match the expansion of an arbitrary theory of gravity at higher orders in hab- 
However, this is not needed to reproduce the first line, i.e., the universal parts up to 0(/i^). 
To wit, the two lines of (B.2) match if we choose 


bi 



1 

4(d - 2) 



r(d/ 2 ) 4-^ \ 
7^+2 pd-l I • 

^AdS / 


(B.4) 


Note that this choice is schematically of the Gauss-Bonnet form. We can interpret this as 
the statement that a Gauss-Bonnet term is all that is needed to reproduce the linearized 
graviton action with arbitrary graviton normalization. This is a generalization of the results 
of the main text: in GFT language, a quantity which is sensitive to both and Cj, may 
undeniably detect the presence of higher curvature terms, but it does not distinguish between 
Gauss-Bonnet theory and any other higher curvature theory in the same way that quantities 
only knowing about a*^ can always be computed using Einstein gravity. The reason is that 
linearized actions up to 0{h'^) only require a small number of terms in the background field 
expansion. 

As a consistency check, we can make the relation between parameters {aoj ^ 2 ) ^ 2 ; ^ 3 } 

and Cj, manifest by referring to the results of [ 68 ] (see also [ 6 ]). There it was shown that an 
action of the form (B.2) leads to \h = cl 2 + 4((i — 2)63 which is proportional to C^, via (2.33). 
This is in perfect agreement with the above statements. 
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Since (B.2) is manifestly written as both a perturbative expansion in the metric fluctuation 
and as a functional of curvatures, we can easily compute the associated perturbative expansion 
of Wald entropy density by taking derivatives with respect to Rated- Note that the way the 
expansion is written makes it explicit that any variation with respect to Rated reduces the 
order in hat by precisely one and the number of derivatives by two. The unperturbed Wald 
entropy is hence only sensitive to and is readily checked to agree with (2.5): 


^Waldi^'n^in-) — 27r 




Vi 


5C 


(I) 

AR 


^ Rated 


^ab^cd 


J AdS 


= W2) a; , (B,5) 


'•AdS 


where the subscript ‘AdS’ refers to evaluation on the background gat- 

Similarly, the perturbation of Wald entropy is now computable via the two contributions 


5Sy/aldi^£n 


= -2tt- 


1 


2t 


id—1 




^ Rated 


R-atR-ed 


(B.6) 


J0(/i) 


where ‘0(/i)’ picks out the term proportional to hat- We find the following contribution from 


C 


(4 . 

AR- 


2tt 


^^Wald, {l)i£min) — 27r 




nd-l 

f^AdS 


(B.7) 


which could also be obtained by simply perturbing the volume form in (B.5). Similarly, one 
can compute the higher curvature contributions from at 0(/i) as follows: 


^^Wald, {II) i^min) 

1 

= -27r' 


2t 


d-i 


Vf 26i - AARef + /^f>ated 

1 


^ab^ed 


= 47r6i4ds-i/ 

^min 

= 47r6ifidsi / + 

S-min 


i? + 2 naeUb^ - ^ Uatriad 


0{h) 


+ Oih^) 


(B.8) 


where is the linearized intrinsic Ricci scalar on Emin in units of ^AdS (i-6., (B.8) is essentially 
the first order perturbation of Wald entropy in pure Gauss-Bonnet theory). In the last step 
we used the fact that extrinsic curvatures vanish on Emin, so the projected Ricci scalar in the 
penultimate line coincides with the intrinsic one. We can further simplify the expression using 

a>7^ = + v„(...) = r^ha^p + V„(...), (B.9) 

^AdS 

Our conventions are as follows: the embedded version of the intrinsic metric pap is Pab = Qab+tatb — UaUb 
where 4 and n“ are normal to the minimal surface such that — —1, rE = 1, t^Ua = 0 and the binormal 

-enK 
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where is the unperturbed intrinsic Ricci scalar and Vq is the intrinsic covariant derivative 
compatible with s-iicl hence only yields boundary terms which we are not keeping track 
of.^* Hence (B.8) becomes 


<55- 


Wald. 


{II) — 27r I A/j 


r(d/2 )e 


TT' 


d/2 


d-1 N 
'^AdS y 




nd—1 


(B.IO) 


Finally adding the two contributions to perturbed Wald entropy, (B.7) and (B.IO), we conclude 
that (B.6) reads 


S S\Y ald{£min') — 27r 



le Vl 

^min * _ ^ 

pd-l 

^AdS 


(B.ll) 


C First law of entanglement 

In the main text we have focused on perturbations of entanglement entropy due to scalar 
deformations. Let us now drop such particular assumptions and consider a slightly more 
general setup which has been argued to serve as a quantitative derivation of bulk dynamics 
from entanglement. For theories in which the Ryu-Takayanagi proposal is valid, a first law of 
entanglement entropy has been shown to be equivalent to linearized bulk equations of motion 
[4-6, 25, 27, 28]. 

The entanglement first law 

6See{A) = Tt{5pa Ha) = 5{Ha) (C.l) 

is purely a CFT statement which says that any change of See{A) due to some small pertur¬ 
bation 5pA of the density matrix is given by the change of modular energy associated with 
the spatial sphere A. (This statement can easily be derived from the formalism of §2.2.) For 
CFTs with a semiclassical gravity dual, this statement has a natural analog in gravity: the left 
hand side of (C.l) corresponds to a perturbation of Wald entropy as in (2.34). Similarly the 
right hand side of (C.l) translates to a change of gravitational energy Egrav{A) of the Rindler 
patch of AdSrf+i associated with the causal diamond of A. That is, the CFT entanglement 
first law (C.l) is equivalent to the gravitational first law 

^S\yfiidi^Smin) — ^Hg^avi^A^ . (C.2) 

The gravitational energy Egrav{A) is defined as an integral of the holographic stress tensor 
over a spatial bulk surface of bulk codimension two, which lies inside the Rindler patch 
of A and ends on dA. Without loss of generality, we push this bulk surface to the boundary, 
so the gravitational energy can be computed as a boundary integral over A'. 



Boundary terms are important to get a finite result as in (2.16). Since we are only interested in identifying 
the fate of the dimensionless parameters, we are being implicit about these subtleties and assume that there 
always exist suitable Gibbons-Hawking-York boundary terms. 
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where is the volume form on the [d — l)-dimensional surface A and C,'^ is the conformal 
Killing vector generating the Rindler flow inside the causal development of A. In holographic 
CFTs the CFT stress tensor of (2.13) agrees with the holographic stress tensor 
One can easily verify that the gravitational energy (C.3) then agrees with the definition of 
modular energy in the CFT, (2.12). 

The holographic stress tensor is determined by the dominant asymptotics of the graviton 
fluctuation on AdS^+i. Let us work in Poincare coordinates 

p2 

ds^ = {dz^ + 77 ^ 1 / dx^dx'^) + hab dx^'dx^ . (C.4) 

Then the relevant graviton falloff near the boundary is given by the following piece of the 
Fefferman-Graham expansion: 

habix^, z) dx°‘dx^ = £^^3 (x^) dx^dx^ -|- ... . (C.5) 

The power of 2 ; is chosen such that the induced boundary stress tensor is finite (i.e., 

2 :-independent). It was shown in [6, 68] that for a wide class of theories the holographic stress 
tensor associated with this bulk perturbation evaluates to 

d C 

= -jf , (C.6) 

where fd is the number defined in (2.33). Note that this form of is consistent with our 

considerations in §2.3.3. Indeed, dSwaldi^min) in the present context can be obtained from 
the gravitational energy integral over the combination (C.6): 

6Swald{£rmn) = ^ ' (^.7) 

This has the same structure as (2.36), i.e., an integral over a quantity (C^, or times a metric 
perturbation) which satisfies a universal equation of motion. 

The linearized equations of motion can be derived from the first law statements above. 
Indeed, assuming the Ryu-Takayanagi conjecture, it was shown in [6, 68] that (C.2) for small 
spherical regions A is equivalent to linearized gravitational equations of motion for any higher 
derivative Lagrangian Cgrav = CgraviOab, Rabcd,'^eRabcd, ■ ■ ■)■ As argued in §2.3, on an AdS 
background these equations are the same for any theory of gravity with second order equations 
of motion, up to overall normalizations which encode the higher curvature contributions. De¬ 
pending on the precise setup, even these normalizations may be absorbed into renormalizations 
of dimensionless quantities. In particular the normalization Xh inside Cj, always multiplies 
(■^Ads/^p)^”^) c.f., (2.33). Unless ^Ads/^P has already been fixed (e.g., by a condition such as 
(2.7)), this ratio can serve to account for any Xh dependence. 
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D Second order gravitational perturbation theory 


In this appendix we list some relevant formulae for gravitational perturbation theory up to 
second order in the graviton. We take the ansatz 


3ab — 9ab hab 


(D.l) 


and are interested in perturbations up to 0(/i^). The background solution gab is AdS^+i 
with curvature — l/^AdS> following simplifications apply for the background curvature 

tensors: 


^abcd — n2 {dacdbd 9ad9bc) j Rab ^2 3ab ; 
^AdS 


^AdS 


R = - 


d{d+l) 

^AdS 


(D.2) 


Let us evaluate and define various quantities, where superscript numbers refer to the order 
in perturbation theory: 

gab = ga^ - ha^ + ha^hl + 0{h^) , 

Kb = Kb + '"'Kb-K «r^, + o(h3), 

oa _ 3a _j_ (i)E)a 

^ bed — ^ bed + ^ bed 

- 'hi ^^R%ad + r%9 


de ^ eb 


+ 0(/l3), 


Rab = Rab + ^^^Rab " Ul ^^^R\cb + r%9 + (D{h^) 


R = R + ^^'>R+-{ ha^hl Rab - ha^ ^^^Rab 

'hi ^^^R\cb + r%g 


gab 


be 


+ (D{h^ 


(D.3) 


In these expressions the linearized objects are given in terms of hab by 

1 


[Vahbd + ^bhad - 
^^^R\cd = V[,Vrf]h^ + VySbhl^ - 






(D.4) 


^^^Rab = VcV(ah^) - -Uhab - -j^a^bh , 

^^^R = ^^^Rab r’^ - R'^^ hab = VaVbha’^ -Oh- Ra^ hab , 
where h = hi = ga^hab- 

For the expansion of actions it is also useful to note the following identity: 




(D.5) 
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